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Abstract. We consider a system of two coupled ordinary differential equa- 
tions which appears as an envelope equation in Bose-Einstein Condensation. 
This system can be viewed as a nonlinear extension of the celebrated model 
introduced by Landau and Zener. We show how the nonlinear system may 
appear from different physical models. We focus our attention on the large 
time behavior of the solution. We show the existence of a nonlinear scatter- 
ing operator, which is reminiscent of long range scattering for the nonlinear 
Schrodingcr equation, and which can be compared with its linear counterpart. 



1. Introduction 

1.1. Physical motivation. In the 30's (see jTO] and [23]), Landau and Zener have 
studied independently the system 

(1.1) — id s u = V(s, z)u ;u(0) = uq, 

where u = t {u\, U2) € C 2 , and the potential V is given by 



V(s,z) = 



s z 
z —s 



This system is the prototype of eigenvalue crossings: the eigenvalues of the matrix 
V(s, z) are V s 2 + z 2 and — Vs 2 + z 2 and they cross when s = and z = 0. More 
recently, a nonlinear version of this system has been introduced in [3] , of the form 

(1.2) 4 ( Ul ) = m ( Ul 

y ' dt\u 2 J y,J \u 2 

with the Hamiltonian given by 

(1.3) ff( 7) = (^) + *(W 2 -H 2 ) * 

where 7(f) = at denotes the level separation, z is the coupling constant between 
the two levels, and S is a nonlinear parameter describing the interaction. This 
nonlinear two-level model can be used to understand Landau-Zener tunneling of 
a Bosc-Einstein condensate between Bloch bands in an optical lattice: this has 
been achieved in e.g. [7]. In Section [2] we present a derivation of this model from 
the nonlinear Schrodinger equation in a rotating frame and a periodic potential. 
Physical properties of the above system have been investigated in e.g. [TBJ [TH1 HZ1 
[20] . We also show how this model can arise under the influence of a double- well 
potential, as in [17j (see Section [2J for a rapid presentation, and the appendix for a 
rigorous proof). The goal of this paper is to recast this model in a mathematical 
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background, and to study some of its properties, in particular as far as the large 
time regime is concerned. 

1.2. Mathematical setting. In the linear case (|1.1|) . classification of crossings, as 
performed in [5] and [B] in the one hand, and |11] on the other hand, produces the 
model problem 

(1.4) -id.u=L^ G ^ju; u(0) = u Q eH 2 , 

where H is a Hilbert space and G an operator on W. Equation (|1.2|) corresponds 
to W = C and G = z G E, as in [115] and [53] (see also [IT] and [5]). Other choices 
of the pair {%, G) are relevant: 

(1) In :6; and [14 , G = z x + iz 2 and U = C. 

(2) In [6] and [9], G = d z - z and H = L 2 (R). 

(3) In [H], [13] and [8], G = q{z) is a quaternion, z G R 4 and "H = C 2 . 

(4) In [8] and [11] , G is a semiclassical pseudodifferential operator on the space 
U = L 2 {R k ),'k G N*. 

For this reason, we will focus on the following abstract problem 

(1.5) — id s u — {^q* u + SF(u)u; u(Q) = u e Ti 2 , 

which is the nonlinear counterpart of (|1.4j) . The nonlinearity F : C 2 —> R 2x2 is of 
the form 

F(u) = diag(F 1 (u),F 2 (u)) 
where the functions Fi,F 2 are gauge invariant, and more precisely, 
(1-6) F J (u) = f J (\\u 1 \\ 2 n ,\\u 2 \\ 2 i ), 

with fj e C 2 (R 2 ; R) and V 2 /j bounded. In particular, the system ([L"2 ]H |1.3 p enters 
into this framework. By the change of variable s = t^/a, we are left with (|1.5|) with 
G = z/^ and Fj(u) = {~iyda-^ 2 (\u 2 \ 2 - \ Ul \ 2 ) for j £ {1,2}. 
We make the following assumption on the operator G: 

(1.7) X (GG*)G = G X (G*G) and x(G*G)G* = G* X (GG*). 

This assumption is satisfied in the first three examples above. The situation is more 
complicated in the fourth one. 

We also assume that the domains of GG* and of G*G - V(GG*) and V{G*G) 
- are dense subsets of H. Note that the domain of the evolution operators q itGG 
and e lrG G are well-defined for any r e R with domain %. 



(IM* + IM*) = o. 



Lemma 1.1. LetUQ G Ti 2 . Under the above assumptions, (|1.5p has a unique, global 
solution u G C(R;'H 2 ). It satisfies the following conservation law 

d 
ds 

Proof. Denote by U(s 2 ,si) the linear operator which maps w lln (si) to u lln (s2), 
where u lln solves the linear equation 

-•".«"' = ( G S . G >' 
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It satisfies U(s 7 s) = Id, U (s , t)U "(t ■, a) — U(s,a), and is unitary on H 2 . By 
Duhamel's principle, ()1.5j) becomes 



u(s) = f/(s, 0)u + »5 / U(s, a) (F(u)u) {a)da. 
Jo 

Local existence follows by a standard fixed point argument (Cauchy-Lipschitz), 
with a local existence time which depends only on HmoII-hxh- We then note the 
conservation law announced in the lemma, which follows from the fact that the 
functions Fj are real-valued. This implies global existence, and the lemma. □ 

We normalize the data so that 
(1.8) IM* + IM* = 1, VsGM. 

In this article, we prove a scattering result for initial data which are in the 
range of lv(o) 2 ^R f° r some R > 0. More precisely, we introduce a cut-off operator 
depending on G: let 9 be a smooth cut-off function , 9 £ Cq°(R) with ^ 9 ^ 1, 
9{u) = for \u\ > 1 and 9{u) = 1 for |u| < 1/2. Then, for R > we set 

9 fl = diag ( e (^-),9 fG * G 



\ r? J> \ R 2 

Because of (JTTTJ) , the operator 9^ commutes with V(s) for all s £ M and QrV(s) is 
a bounded operator on "H x H with norm Vs 2 + R 2 ■ Besides, a simple computation 
shows that m_r(s) = Qrii(s) satisfies id s UR = V{s)ur + F(u)ur with ur(0) = 
QrUo = wo- Therefore, we have the following result. 

Lemma 1.2. Suppose uq — QrUq for some R > 0, then for all s £ M, £/ie solution 
of (jl.5l) satisfies u(s) = Oru(s). 

Typically, in the physical examples presented in Section [2j the assumption uo = 
Qruq consists in saying that some physical parameter (whose value is fixed in 
practise) belongs to a bounded set. 

1.3. Scattering. We introduce the phase function cp given by 

s 2 A 
p(s,A) = — + 2 ln l s l- 

We can describe the large time asymptotics of u. 

Theorem 1.3. Assume that u$ — QrUq for some R > 0. Then, there exist a — 
(0(1,0:2) £ H 2 , cu = (wijWq) G H 2 , such that: 
I. As s goes to —00. 



u 2 (s) = e^W-WOo, + o {R 2 



2. .As s goes to +00. 

U 2 ( a ) = e M*5i(«-)»-< V (-.GfGf) W2 + (i? 2 S -!) 
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The above result is reminiscent of long range scattering in nonlinear Schrodinger 
equations, as described first in [21] : nonlinear effects are present both in the fact 
that the amplitude of the functions Uj undergoes a nonlinear influence (this is what 
happens in nonlinear scattering in general), and in the fact that oscillations are 
different from those of the linear case (a typical feature of long range scattering), 
since the linear phase ip is not enough to describe large time oscillations. An 
important difference though is that (|1.5j) is not a dispersive equation, as can be 
seen from (11.81) . 

Remark 1.4. When 6 = and G is scalar (G = z), this result goes back to the 30's 
with the proofs of Landau and Zener [T5] and [23] . The original proof is based on 
the use of special functions; more recently, the proof of [TU] relies on the analysis 
of oscillatory integrals. When G is operator- valued, the theorem is proved in the 
linear case (S = 0) in [T2J Proposition 7]. We point out that there is a slight 
difference with the present framework, due to nonlinear effects. In the linear case, 
one associates with uq scattering states such that the asymptotics hold true for 
&ru(s). Here, the scattering states depend on R in a non trivial way. 

Conversely, wave operators are well-defined, as stated in the following result. 

Proposition 1.5. Let a = (01,02) £ Ti 2 such that a = Qua for some R > 0. 
There exists C > 0, such that for all e > 0, there exists a solution u e £ C(M.;H 2 ) 
to (fi~5l) with 



lim sup 



H 



< Ce. 



Remark 1.6. A similar result holds in +oo. 



Therefore, there exists K, C % a dense subset of 1-L and an operator Ss with 
domain K, which is a nonlinear scattering operator mapping the scattering state 
at -co, a £ K,, to its counterpart uj at +oo given by Theorem 11.31 the scattering 
operator is defined by ui = Sga. 

Besides, if the unit ball of T-L is a finite dimension vector space (typically, T~L = C 
in the case of (|1.2I) ). one can extract from u £ a converging sequence with limit u 
which solves (( 1 . 5|) and admits the scattering state ( 0:1,0:2), and the operator S$ is 
defined on H (fC = %). 

Our final result concerning the large time behavior of solutions to (|1.5p consists 
in describing the effect of the nonlinearity in Sg by comparing this operator with 
its linear counterpart. We denote by u lln the solution of (|1.4[) . and we associate 
with uo the linear scattering states 

~ lin /~lin „lin\ , .lin /, ,lin , ,lin\ 

a = (oi , a 2 ) and u = (w x , ui 2 ) . 
According to Proposition 7 in |12j . the linear scattering operator is given by 



(1.9) 
with 



lUn _ / a(GG*) -b(GG*)G 
' ' \b(G*G)G* a(G*G) 



«(A)=e-*, HA) = ^2-«/V^r(l +i ^) sinh(^ 

\2 



and a(\y + X\b{\)\ 2 = 1. 
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When G — z, the coefficient 

(1.10) T{z) = a(z 2 ) 2 = e"^ 

is the celebrated Landau-Zener transition coefficient which describes the ratio | 2 /l a i | 2 
of the energy which remains on the first component (when a!, 111 = 0). As we shall 
see in the next result, the Landau-Zener transition probability remains relevant in 
the nonlinear regime and for small 5. Let us define, for j = 1,2, 

(1.11) At = / (F 3 (u lin (r))-F 3 (^))dr 7 

Jo 

(1.12) AT = / (Fj (u lin (r)) - Fj (a 1 ™)) dr, 

Jo 

and A* = diag(Af , A*). 

Theorem 1.7. Let S > 0. TTie nonlinear scattering operator S$ and the linear 
scattering operator S hn are related by the formula 

S s (a) = e * 5A+ («)s lin (V l5A ~^a) . 

In addition, we have the following asymptotic behavior as S —¥ 0. There exists a 
constant C > such that for all R > and initial data uq with Qruq = uq, we 
have 

\\S s -S lin -iS (A+S lin -S lin A-)\\ c(Hxn) ^CRS 2 . 

As expected, as 5 — > 0, Ss behaves at leading order like the linear scattering 
operator. Nonlinear effects show up in the 0(§) corrector (the nonlinearity is 
present in the definition of A~ and A + ), and are described rather explicitly in 
terms of linear scattering and of the nonlinearity F. We note also that if A^ ^ k\ 
{Fx 7^ F 2 ), then S hn and A+ do not commute. 

This paper is organized as follows. In Section[21 we sketch the derivation of mod- 
els of the form (11.5[) from cubic nonlinear Schrbdinger equations used to describe 
Bose-Einstein Condensation. In Section [3j we set up some technical tools needed 
for the large time study of (|1.5j) . and we prove Lemma ll.2l Theorem ll.3l is proved 
in Section |H Proposition [L5] in Section [5j and Theorem 11.71 in Section HI Finally, 
in Appendix [S] we go back to the derivation of (| 1 . 5|) from physical models, and 
establish rigorously that (jl.5l) can be interpreted as an envelope equation in the 
semi-classical limit. 



2. Formal derivation of the model 

We rapidly describe some cases where (jl.5p appears as an approximation to 
describe the motion of a Bose-Einstein condensate. 

2.1. Condensate in an accelerated optical lattice. As proposed in j3], and 
considered in e.g. |16j . the motion of a Bose-Einstein condensate in an acceler- 
ated ID optical lattice is described by the equation 

ih-^- = — ~ih— -ujt) Tp + Vo cos (2k L x) t/> + -|^| V, 
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COS 



2 Akx 



where m is the atomic mass, kh is the optical lattice wave number, Vq is the strength 
of the periodic potential depth, u) is the inertial force, and a s is the scattering length. 
The adimensionalized form of this equation reads 

dib if d \ 2 

(2.1) i-jZ = - f -i—-at\ ip + vcos(x)ip + e\i()\ 2 ip, 

with e = —1 or +1 according to the chemical element considered. The approach 
in [3] consists in injecting the ansatz 

(2.2) i/j(t, x) = a{t)c lkx + b(t)c l{k - 1)x 

into (|2.1j) . with k = Ul = 1/2, corresponding to the Brillouin zone edge (this 
approximation amounts to considering that only the ground state and the first 
excited state are populated, see 16 ). We compute 

id t rp = iae ikx + ibe^ k -^ x , 
(-id x - at) 2 $ = (k - at) 2 ae ikx + (k - 1 - at) 2 be l{k ~ 1)x , 

<x)ip = ^ (ae 4(fe+1)a; + be lkx + ae l{k ~^ x + be l{k ~ 2)x ^j , 

= (| a | 2 + \b\ 2 ) (ae lkx + be l(k ^ x ^j + |a| 2 6e l(fe - 1)a: + a|6| 2 e 

+ a 2£ e i(fc+l)x +m 2 e i{k-2) X _ 

Leaving out the new harmonics (the last two exponentials) generated both by band 
interaction and nonlinear effects, and identifying the coefficients of e lkx and e^ k ~ l ' x , 
we come up with: 

id t a = - (k - at) 2 a+^b + e (\a\ 2 + 2\b\ 2 ) a, 
id t b =-(k-l- at) 2 b+^a + e (2\a\ 2 + \b\ 2 ) b. 

We notice the identity d t (|a| 2 + |6| 2 ) = 0, so we can write \a\ 2 + \b\ 2 — m 2 > 
(in [3], mo = 1). Now recalling the numerical value k = 1/2 and expanding the 
squares, we have 

1 at (at) 2 v, , 2 ,,. 2 x 
id t a = -a- —a + ^r^a + -b + e (m£ + |6| 2 ) a, 
o Z Z Z 

ldth = 8 T 2 + 2° + 6 ( a + |a| ) 
The above system becomes 

* = { E ° + + 5 (T i) + 1 Di - g w- |o|J > (i 

with 

1 9 

= g + mo- 
using finally the gauge transform 

Ul ^\ = e -iE t-ia 2 t 3 /6 (a 
u 2 
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we end up with 



Ui \ 1 ( —at v \ ( u\\ x e TA . _ n |2 , l2 \ ( u i 



which is of the form (|1.2[) - (|1.3|) or (|1.5|) via the change of variable s = i^/f ; we 
then have G = (2a)~ 1 / 2 v and F(u) = e(2a)- 1 / 2 Diag(|u 2 | 2 , K| 2 )- Note that the 
only approximation that we have made in this computation consists in neglecting 
the other harmonics than q ±ix / 2 ^ and no linearization was performed, unlike in the 
computations of [3]. 

Theorem 11.31 gives an asymptotic for the profiles a and b of the ansatz (|2.2|) as t 
goes to ±oo. Theorem 11.71 gives an information on the profiles (a+, 6+) for t ~ +oo 
in terms of (a_,6_), those for t ~ — oo. For example if for t ~ — oo, we have 
(a_,6_) = (a, 0), then the profiles for t ~ +cxi are related via the Landau-Zcncr 
transition coefficient (|1.10j) for z = (2a)~ 1 / 2 v: at leading order in <5, they satisfy 



|a| 2 and |6+| 2 = (l-e^) \a\ 2 . 



( 2 - 5 ) w = w = X;,? . / 3 >°- 



2.2. Condensate in a double- well potential. As suggested in [3], and further 
developed in |17j . nonlinear Landau-Zener tunneling may be realized in a double- 
well potential. Consider 

(2.3) ^ + 10 = W(t)a# + e |^|Y 
The potential W is of the form 

(2.4) W(t,x) = V s {x)+KtV a (x), 

where V s is symmetric and V a is antisymmetric. In |17j . we find the explicit ex- 
pression 

Vo , _ Fosinh^x) 

cosh 08a:)' MXJ cosh 2 (^) 

Note that V s is a double-well potential. The main point to be aware of is that the 
lowest two eigenvalues A + < A_ of the Hamiltonian — \d 2 + V s are non-degenerate, 
with associated eigenfunctions ip± (see Appendix A for details) . The two exponen- 
tial functions e ±lx / 2 of the above model are then replaced by the so-called single- well 
states 

(2-6) (f L = -7= (<f+ ~ <P-) , <fR = -y=(<p+ + <p-). 

We note that this approach can be generalized to a multidimensional framework, in 
the spirit of [23] (see Appendix |A"|) . We sketch the computation in the simplest ID 
case though, to emphasize the differences with the optical lattice case. We shall 
use mostly two properties of ip^ and tp R , described more precisely in Appendix [Al 

• Vr{~x) = Vl(x). 

• The product tpLfR IS negligible, because ipL and tpR are localized at the 
two distinct minima of V s . 

Seek ijj of the form 

ij)(t, x) = a L (t)ip L (x) + a R (t)ip R (x). 

Denoting 

A_ + A+ A_ - A+ 



fl = 



LU = 



s 
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we compute: 

id t ip = ia L ip L + ia R ip R , 

~\® 2 j^ + = aL i^VL - uipb) + a R (iltfR - lolpl) 
= (ila L - wa R )tp L + (tt R - uja L )tp R , 
V a tp = a L V a ip L + a R V a ip R , 
\ip\ 2 i> = (|a L | 2 |v3 L | 2 + la^l 2 !^! 2 + 2Re (a L a R lp L ip R )) (a L (p L + a R ip R ) . 
By integrating in space and neglecting the product Pl^Pr, we get: 

{idtaL = SIcll - u;a R + Kt^/LaL + eL\a,L\ 2 aL, 
id t a R = Vla R - oja L + nt-f R a R + e R \a R \ 2 a R , 

with 

C A „ f T/ ,„2 f 4 



1L = / V a <p L , e L = e <p L , j R = V a ip R , e R = e ip R . 

JR JR JR JR 

By symmetry, ji — —"f Rl £l = e R =: 5, so if we set a = kjl, we come up with: 
J id t a L = Qa L - usa R + ata L + S\a L \ 2 a Ll 
[ id t a R = fla R - cuaL - ata R + 6\a R \ 2 a R . 

Using the gauge transform 

ui\ ( a_Le* n * 



u 2 J \a R e" nt 

we find: 

.„ fui\ fat —(jj\(ui\ lX (\ui\* \ /hi' 

^UJ = U> -<*){<») +6 { o M 2 JU 

which is of the form (|1.5|) via the change of variables s — y/at with Fj(u) — 
a^ x / 2 \uj\ 2 and G = -ua -1 ' 2 . Therefore, Theorems 11.31 and 11.71 yield similar 
results than in the preceding subsection with a Landau-Zener coefficient e~ n ~ . 

3. Technical preliminaries 

In this paragraph, we introduce several operators related with G that will be 
useful in the following. We set 

(3.1) J := (l and K := V(0) = % 



V(s) = sJ + K. 



so that V(s) writes 
We observe that 

V(s) 2 =A( S ) 2 , 
where A(s) is the diagonal operator 

A(s) = diag s 2 + GG*, \/ s 2 + G*G*) . 

For this reason, A(s) appears like a diagonalisation of V(s), all the more that if we 
set 

n ± (.s) = i(ld±A( s )-V( s )), 
then we have the following properties. 
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(1) n±(s)F(s) = ^(s)n±( s ) = ±A( s )n±( s ) = ±n ± (s)A(s). 

(2) n+(s) + n-(s) = w. 

(3) (n±( s ))* = n±( s ). 

(4) n±(s)nT( s ) = o and (n±( s )) 2 = n±( s ). 

The properties (2)-(4) show that n ± (s) are orthogonal projectors, and the prop- 
erty (1) will play the role of a diagonalisation of the operator V(s). The fact that 
II ± (s) and V(s) commute with A(s) is more general. In fact, V(s) and n ± (s) are in 
the subset A of C(H 2 ) defined by: A £ A if and only if there exist smooth functions 
a, b, c and d such that A = A(a, b, c, d) with 

a{GG*) b(GG*)G\ 
c{G*G)G* d(G*G) ) ' 

A simple calculation shows that, because of the commutation property (|1.7j) . oper- 
ators of A commutes with A(s): 

\/A G A, V.s G M, AA(s) = A(s)A. 

Besides A is an algebra, as shown by the following lemma which stems from straight- 
forward computations 

Lemma 3.1. Let a,b,c,d,a' ,b' ,c' ,d' G C°°(M). We have 

A(a, b, c, d)* = A(a, c, &, d), 
A(a, b, c, d)A{a', b', c', d') = A(a", b", c" , d"), 

with 

a" (A) = a(A)a'(A) + A6(A)c(A), 6"(A) = a(A)6'(A) + b(X)d'(X), 
d"(X) = d(X)d'(X) + Xc(X)b(X), c"(A) = a'(A)c(A) + d(X)c'{X). 
Operators of A will be called diagonal if 

A = U+(s)AU + (s) + U-(s)AU-(s), 

and antidiagonal if 

A = U + (s)AU'(s) +r(s)An + (s). 

In particular, V^s), II + (s) and II - (s) are diagonal operators of A; on the other 
hand, the operators <9 s II + (s) and d s Il~(s) are antidiagonal elements of A. Indeed, 
the relation 

d s U+ = <9 S ((II + ) 2 ) = U+d s U+ + d s U+U+ 

implies that n*^!!" 1 "!!* = (and similarly for II~ since d s II~ = — d s Il + ). 

Antidiagonal operators have nice properties that we shall use later. Typically, 
they can be written as commutators with V(s): if C(s) = IT ± (s)C(s)n =F (s), then 
we have 

C(s) = ±[B(s),V(s)], 

with 

(3.3) B(s) = ±A( S )^C( S ), 

which also belongs to A. Because of this property, we have the following lemma. 



(3.2) 



A(a, b, c, d) 
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Lemma 3.2. Let u be a solution of (11. 4|) . C(s) € A with C(s) = II + (s)C(s)II (s), 
and let B(s) be associated with C(s) as in (|3.3|) . then 

(C(s)u(s) , u{s)) U 2 = ~ (-B(s)u(s) , u(s)) K2 
(3.4) + * (d s B{s)u{s) , u(a)) wa - ([B(«), F(u)]u(s) , u(a)) w „ . 

Proof. We write 

(C(a)«(«) , = ([£(*) , , u(s))„ a 



B(s)i9 s ii(s) , u(a) ) - Ib(s)u(s) , -<9 s u(s) 
H 2 \ 

-([B(s), F(u)]u(s), u(s)) n2 
ld_ 
i ds 



7— ( B ( s ) u ( s ) . u ( s ))n 2 



1 

and the lemma follows. □ 



4. Existence of scattering states 

The proof of the existence of scattering states (Theorem ll.3j) consists in three 
steps. 

(1) We set 

u±(s) = TJ±(s)u{a), 

and we first prove the existence of a limit to ||u ± (s)||^ 2 as s goes to ±00. 

(2) Then, we deduce that for j 6 {1,2}, the functions Fj(u) have a limit that 
we denote by Fj(uj) (resp. Fj(a)) as s goes to +00 (resp. —00). Besides, 
we prove that the functions 

r+OO l-S 

shW (Fj(u(T)-Fj((jj))dT and / (Fj (u(r) - Fj (a)) dr 

J s J — 00 

are well defined and we study their behavior as s goes to +00 or —00, 
respectively. 

(3) Finally, we deduce the existence of scattering states from the scattering 
result for the linear equation and by using Duhamel formula to treat the 
nonlinearity. 

If one expects the asymptotics behavior of u to be described as in Theorem ll.3[ then 
the first two steps consist in deriving the asymptotic phase, while the goal of the 
final step is to describe the amplitude. This strategy is similar to the one employed 
in e.g. |15j to study the long range nonlinear scattering for the one dimensional 
cubic Schrodinger equation and for the Hartree equation. 

First step. The first step of the proof relies on the following proposition. 
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Proposition 4.1. Assume that uq — Qruq for some R > 0. Then, there exist 
il ± , A ± ^ such that 

hH*)\&xH-&=°(Jji) ass^+cx., 
h ± (s)\\ 2 nxH -A ± = o(^j as S ^-oo. 

Proof. We consider the limit s — > +oo for the + mode. The limit s — > — oo and the 
case of the — mode can be treated similarly. We are going to prove that ||u + (s)||| ixW 
is a Cauchy sequence as s — > +oo. In that purpose, we write the equation satisfied 
by u + 

-d s u+ = A(s)u + + -d s U+{s)u + F(u)u+. 

i i 

Therefore, for < t < s, 

(4.1) \\u + (s)\\ 2 H2 -\\u + (t)\\ 2 nxn 



= 2Re^y (U+(a)d s U+(a)u{a), u{a)) UxU da 



We are going to use properties of the operator II + (s)(? s n + (s) that we gather in 
the next lemma where we denote by ImC the skew adjoint part of the operator C: 
ImC= (C-C*)/2. 

Lemma 4.2. Let C(s) = U + (s)d s U+(s) and B{s) = iA(s)- 1 C*(s). Then C*(s) 
and B(s) are antidiagonal operators of A with C(s) — O (-), B(s) = O (p-) and 
d s B(s) =0 (4) in C{H 2 ). Moreover, 

The proof of this lemma is postponed to the end of the section. Using LemmaGO 
we obtain 

(C(a)u(a) , u(a)) H 2 da 
= - (B(a)u(a) , u(a)) nxn 

+ * It ^ sB ^ U ^ ' u ( cr ))« 2 d° _ J (l B ( <T )>F( u )] u ( <7 ) > u { <J ))w da. 
A short computation gives for s£l, 

Bb{[B(s),F(u)]u(s) , u(s)) w =Bje{[lmB(s),F(u)]u(s) , u(s)) w , 

with 

[lmB(s),F(u)} = ^(F 1 (u)~F 2 (u))A(s)- 3 K. 

Therefore, 

Re([B( S ),F(u)}u(s) , u(s)}„ 2 = ^(u) - F 2 (u)) (A(s)~ 3 Ku,u) n2 . 

Since it(s) = 0_Rit(s) with ||w(s)||%2 = 1, we have ||iftt||^a ^ R. We deduce for all 
s G K 

|Re<[B(«),F(u)]u(«) , u( 8 )) wa | =0 ^ '''' 
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This observation and Lemma 14721 yield: for s > t 3> 1, 

\\u + (s)\\ 2 nxn -\\u + mu* = o(§ 
This implies the convergence of ||m + (s)||-^ x ^ to Q + with 

||m+(s)||^ 2 = n+ + o 

hence the proposition. □ 



Proof of Lemma The proof relies on the computation of C (s) by observing 

d s n+( s ) = i (A^)- 1 j - sA( s )- 3 v(s)) , 

where we have used d s A(s) = sA(s) -1 . Since 

n+(s)a s n+(s) = n + (s)9 s n + (s)n-(s), 

the operators C(s) and B(s) are antidiagonal. Besides, we have 
C(s) = U + (s)d s U+{s)U-(s) = ^A(s)- 1 U + (s)JU-(s) 

= iA( S )-V + A( s )" 1 [y( s ) ) J] - A(s)- 2 V(s)JV(s)). 

In view of 

[V(s), J] = [K, J] =2 (q+ ~ G J , 7(s)JV(s) = s 2 J + 2sK + di&g{-GG*,G*G), 



we obtain 



1 w s-3 ^ ; „„^~* „ M| liM-j/0 -G 



C(s) = -A( S )" J (diag(GG* , G*G) - sK) + -A(s) { f , : () 



hence (|4.2[) . Then, the estimates 

(4.3) IIAW^IUcw") < a -1 , II fl.ACs)- 1 ]!^) < 1 
||a s 2 A( s )||« 2 = ||A( s )- 3 diag(GG*,G*G))||« 2 = O^ 1 ) 

allow to conclude the proof of the lemma. □ 

Second step. We study F(u) by analyzing the norm of u%(s) and U2(s). Note 
that, for \s\ 3> 1, we have 

A (s) -e fl =i|(id + (5 

We deduce asymptotics for the operators n ± (s): for s > 0, in £("H 2 ), 

(4.4) n+(.s)e fl = 9 R £i + Ia^KQr + O ( h '~ 



U-(s)Q R = Q R E 2 - ^A^KQr + O f U ~ 
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and for s < 0, in C{% 2 ), 



Il-(s)Q R = QrEi - Ik{s)- 1 KQ r + o(^ 

where K = V(0) has been defined in (|3.1[) . 

E x = diag(l, 0) and E 2 = diag(0, f ). 

This suggests that as s goes to +oo, ui(s) is related to u + (s) and u 2 (s) with u~(s) 
while when s goes to — oo, Ui(s) is related to u~(s) and U2(s) with u + (s). For this 
reason, we denote by F(cj) and F(a) the matrices 

F(uj) = F(fl + ,fl~) and F(a) = F(A~,A + ). 

We now prove that F(lj) and are the limits of F(u) as s goes to ±oo. More 

precisely, we have the following result. 

Proposition 4.3. If OrUq = uq for some R > 0, then for j £ {1, 2}, 



+oo, 



(Fj (u(t)) - Fjtu)) dr = ( — J as s 

£ (F j (u(r))-F j (a))dr = o(j^j as s -+ -oo. 

Proof. We first transfer the result of Proposition 14.11 on the coordinates of u(s). 
Recall that uq = QrUq yields that for all s € M, u(s) = Qru(s) (by Lemma \1. 2 1) . 
We claim that for all R > 0, as s —> +oo, 



(4.5) 



\u 2 (s)\\ 2 n - fl- = g 2 (s)+0 



s- 
FT 

9 / ' 



where for j £ {1, 2}, \ gj (s)\ = C(i? S - 1 ) and J s + °° gj(a)da = 0(R 2 s~ 1 ), and that, 
as s — > — oo, 

R 2 



IM*)llw-A-=5i(*)+0( ^1, 

IM*)ll«->4 + = ^(*) + o(5)' 

where for j e {1,2}, |<fc(s)| = ©(^Is^ 1 ) and ^°° gj (a)da = 0{R 2 \s\- 1 ). 

Then, one concludes as follows: for r > and j € {1,2}, we write 
F 3 (u(r)) - Fj(u) = V/, (Sl+,Sr) ■ (gi(r),g 2 (T)) 

+ f S7 2 f, (n+ + a 9l (r),n- + ag 2 (r)) [ gi (T),g 2 (r)} 2 (1 ~ a)da + 0(R 2 s- 2 ), 



where the fx, f 2 are the functions of (jl.6l) . Recall that, for j £ {1,2}, V 2 /j are 
bounded, V/, are continuous and Qf ^ 1 from (|1.8|) . Therefore, the fact that 
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\gj(s)\ = OiRs- 1 ) and / s + °° 3j (r) = 0{R 2 s- 1 ) for all j e {1, 2} yields the lemma 

for s — > +oo (one argues similarly in — oo). 

Let us now prove (|4.5|) . Indeed, by Lemma l4.1[ there exists fl + such that as 
s — » +oo, 

||m + (s)||^2 = n+ + o(r s - 2 ). 

By (|4.4|) . we can write 

?2> 



2 



= ||« + (*)||^ +Re(n+(s) , A^)" 1 ^^)}^ + O (^j 

= h + (s)\\ 2 n 2 +Rc(KA(s)- 1 Il+{ S )e R u( S ) , u(s)) n2 +o(^ a 
In view of (|4.4p and of the relation 

we have 

n+(.s)e fl ^Sxe^ + o^s- 1 ). 

This implies that we can write 

k A( s )- 1 n+( s )e fl = iTA(s)- 1 E 1 e R + o ( % 



= E 2 KA( s y 1 E 1 e R + o (jjpj 
= A( s )- 1 n-( s )^e fi n+( s ) + o (% 



where we have used E 2 K = KE\, A(s) = £'iA(s) 1 and the commutation 
properties of V(s) and II (s) with A(s) _1 and 0_r. We set 

We obtain 

(4.6) 9l (s) = OiRs- 1 ) ^d\\u 1 ( S )\\ 2 n = \\u+( S )f n2 +g 1 (s) + o(^ 

In order to prove the integrability properties of g\, we study the operator 

C(s) = A(s)- 1 U-(s)Ke R U+(s). 

Lemma 4.4. The operators C {s) = A(s) _1 II- (s)KQ R U + (s) andB(s) = |A(s) -1 C70) 
are antidiagonal operators of A which satisfy 

(4.7) ImC(s) = ^A(s)- 2 e R ^ ~ G 
and C(s) = O (f ) , B{s) = O (J- ) , 0.5(a) = O *n £(H 2 ). 



We postpone the proof of this lemma. By Lemma 13.21 we write 
1 d 



- Re \j— (B(s)u( S ) , u(s))) + Re (i (d s B(s)u(s) , «(«)) wa ) 
-Re([£(s),.F(u)]u(«) , u(s)J) 2 +0' " 
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In view of 
we obtain 



[lmB(s),F(u)] = -(F 1 (u)-F 2 (u))A(s)- 3 K, 



9l (s) = Re [~ (B(s)u(s) , u(s))) + Re (i (d s B(s)u(s) , «(*)) wa ) 

+ |(Fi(u) - F 2 ( u ))Re(A( s )- 3 ifu(,s) , u(s)) n2 + O (^j . 

By ((46]and using that sA(s)~ 3 < s~ 2 , ||iCu(s)||^a < i? and \\u(s)\\ H 2 for all s G K, 
we infer (|4.5|) with the properties stated for gi(s). The other assertions of the claim 
can be proved similarly. □ 



It remains to prove Lemma 14.41 

Proof of Lemma \4-4\ We write 
1 

4 

In view of 

—G 



C(s) = -Hsy^K + A( s )- 1 [ J ftT, V(s)] - A(s)- 2 V(s)KV(s))Q R . 



[K,V(s)]=2s 
and 



G* 



V(s)KV{s) = -s 2 K + 2sdiag(GG*, -G*G) + I qJqq* q 



we obtain 

C(s) = iA(s)- 1 (-2sA(s)- 2 diag(GG*, -G*G) + (1 + s 2 A(s)- 2 )AT 

hence (|4.7I) and the other properties of the Lemma follow from the observation 
||-K"u(s)||^2 = O(R) and from the properties of A(s) stated in (|4.3[) . □ 

Third step. We are now able to prove the existence of scattering states and to 
conclude the proof of Theorem 11.31 



Proof of Theorem \1.3\ Let us focus on the case s — > +oo. Because of Lemma 1473 
we can define v(s) = (i>i(s), V2(s)) with for j = 1,2, 



Vj(s) = Uj(s) exp I — isSFj(uj) — iS J (Fj(u(r)) — Fj(u))di 
We have 

(4.8) -d s v = V(s)v; Wi s=0 = u . 

i 

The linear scattering result of [HI Proposition 7] can then be applied to v. This 
gives the existence of u> = (tDi, 6)2) (independent of R) such that 

Vl (s) = e lvis ' GG ' ] 9(GG* /R 2 )Q! + o(l), 

v 2 {s) = e-^ G ' G '>d(G*G/R 2 )Cj 2 + o(l), 

where 9 is the smooth cut-off function used for defining the operator Or. Besides, 
the rest term o(l) which comes from the o(l) in Lemma 7 in [T2] (or equivalently 
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Lemma 11 in [IT]) has been studied in the proof of Proposition 5.5 of [13] and 
proved to be 0(R 2 s~ 1 ). Therefore, we have 

Vl (s) = e* v ^' GG ^e(GG*/I^)Q 1 + ©(i? 2 ^ 1 ), 

v 2 (s) = e~ l ^ s > G * G) 0{G*G/R 2 )Cj 2 + 0(R 2 S - 1 ), 

Note that v does not need to be truncated by O r since it is associated with a initial 
data which is already truncated by Qr (uq — OrUq). Again by Lemma 14. 3[ we 
have 

exp iS / (Fj(u(T)) - Fj{u:))dT =1 + 0' 



Therefore, 



Ui(s) 
u 2 (s) 



= e 4S<5F2(tj) v 2 (s) = e ^F 2 (u 1 )-t V (s,G*G) uj2 + o{R 2 s - 1 ), 



with 
(4.9) 



u>t = exp [i6j+°° (F x (u(r)) - F x {u))\ 6(GG*/R 2 )u u 

uj 2 = exp [iS J +o ° (F 2 (u(t)) - F 2 (w)j\ 9(G*G/R 2 )u 2 , 
and the existence of scattering states in +oo is proved. We point out that 
IM*)I& = WviWfa, \\u 2 (s)\\ 2 n = \\v 2 (s)\\ 2 n , 

so that we have 

n+ = \\6{GG*/R 2 )u l \\ 2 H = \\uj l \\ 2 n , 
Q- = WO^G/R 2 )^ = Hwall^, 

which justifies the notation F(u>) — F(tt + , f2~). 

A similar argument allows to define the scattering states in — oo. 



□ 



5. Existence of wave operators 

In this section, we prove Proposition 11.51 Let a = (ai,a 2 ) £ Ti 2 such that 
a x £ V{GG*) and a 2 £ V{G*G) with 

ll"i||« + ll"2||« < 1 and \\GG*a 1 \\ H + \\G*Ga 2 \\ n ^ R 2 
for some R > 0. This assumption implies that ||G*ai||% + ||Gq!2|| / h ^ R- 

The strategy is the following. We first introduce an ansatz which solves (II .5|) 
up to an integrable source term. Then, for e > 0, we infer the existence of an 
exact solution which is at distance e > of this ansatz as s — > — oo. We conclude 
by noticing that the ansatz is an 0(1/ \s\) perturbation of the expected reference 
solution. 

Step 1. Construction of the ansatz. We consider a continuous function from 
R to H 2 , u(s) = (ui(s),u 2 (s)), with 



ui(s) 
u 2 {s) 



2(«) 



a 2 



s 

je^s) (fa 



7i 



72 
c 2 



-is 1 /I 



S 2 

«a(a) . 



3 3is 2 /2 



-3is 2 /2 



«l(g) 

S £ 
K 2 (s) 
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The phase operators B\ and 2 correspond to the phases obtained in the previous 
section, and are given by 

9i(s) = y + M*), M*) = lQn\s\)GG* +6sF 1 (a), 

9 2 (s)=-y + Ms), Ms) = -l(Ms\)G*G + 5sF 2 (a). 
Similarly, we dchne 

e 1 (s) = ^- + ±(ln\s\)G*G + 5sF 1 (a), 
9 2 (s) = - S - - ~(ln|s|)GG* + 6sF 2 (a), 

so that we have 

(5.1) e^G* = G*e lf>1 and Ge* = e Wl G, 
e*G = Ge* and G*e i0 ~ 2 = e z0 *G*. 

The coefficients j3j, jj and Kj, j € {1,2} satisfy 

j5 2 = -jG*ax, (3i = --Ga 2 , 

s s 

7i = ^(F 2 (a) - Fiiaftpi, 72 = ^(F 2 (a) - Fi(a))0 2) 
= -|e'*ai (ft/i(a) ^e^^ GG ' ai )^ + d 2 h{a) (e-^^ G ' G a 2 , ft) J , 
«a(») = ^ 2 a 2 (ft/ 2 (a) (e l(ln|s|)GG *c*i, ft) + ft/ 2 (a) (ft, e- l ( ln l s l> G * G a 2 ) J , 

Sa(«) = "^ 2 a 2 (ft /a (a) (ft, e 4 ( ln l s D G * G ai )^ + ft/ 2 (a) ( e -^ ln l s l> G * G a 2 , ft) J , 
with 

^>i = ft - - *i(a)) = i(ln|«|)GG* + SsF 2 {a), 

fa = fa + 5s(.ft(a) - F 2 (a)) = i(ln|s|)GG* + *«(2f 1 (a) - ft(a)), 

and 

fa = 2 + 5a(Fi(a) - F 2 (a)) = --(ln|s|)G*G + 5sF 1 (a), 

ft = ft - 5a(Fi(a) - F 2 (a)) = --(ln|s|)G*G + £s(2F 2 ft) - F 1 (a)). 

We have denoted by difj(a) the real-number 9j/j(||ai||^, ||a 2 ||^)- 
We observe that for j £ {1, 2}, 

(5.2) ||ft|| H = 0(E) and || 7j |ft = (fti?), 

Il«i(*)ll7i + = O(ii) and ||ft^(s)||« + ||ft^(s)||« = (fti? 3 ^ 1 ). 

We begin by studying each term of the equation (|1.5|) : s-derivatives, linear term, 
nonlinear term. Then we project on each family of oscillations generated by the 
operator- valued phases. 
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Derivatives. Using f|5 . 2[) . we get 

= e l(>1 (s + y s GG * + SF^a^j ai + e i§ * (-fa + ^F 2 (a)fa - ^ 

- -e-^Kiis) + -e^nUs) + C(i? 3 s - 2 ), 
s s 

-id s u 2 = e 102 (s - ^-G*G + SF 2 (a)^j a 2 + e^ 1 (fa + 5 -F x (cc)fa + y 

+ -e l ^K 2 {s) - -e" 3! ^ 2 ( S ) + 0(R 3 s- 2 ). 
s s 

Linear part. Using (|5.ip and (|5.2p , we get 

sui + Gu 2 = e tf>1 (sai + - g Gfa^j + e 102 (fa + ~ 7l + Ga 2 

+ - e - i 'r Kl f 8 ) + -e 3i ^Rx(s) + 0(R 2 s- 2 ), 
s s 

-su 2 + G*u x = q 102 (~sa 2 + ^G*fa^j + e* 91 (-fa - \ 2 + G*«i 

_ I e l ^K 2 (s) - - e -^k 2 (s) + 0(R 2 s- 2 ). 
s s 

Nonlinear part. We begin by expanding Fj(u) = /j(||tti||^, ||w2||%) for j e {1,2}. 
We observe 

ll«i||« = l|ai||« + -Re(e i8l a 1 ,e i§2 fa) + 0(R 2 S - 2 ), 
2 



\u 2 \\ 2 H = \\a 2 \\ 2 H + - Re (e l02 a 2 , fa) + 0(R 2 S - 2 ). 
s \ / « 



We deduce 



(5.3) 



Fj(u) = Fj{a) + Re (e i8l ai ,e i82 fa) 

2 



+ ~9 2 / J -(a)Re(e i92 a 2 ,e i9l /3 2 ) + 0(i? 2 S - 2 ). 



e -i9 2e i6li = e is 2 +i5 s (_F 1 ( Q )^F 2 (Q)) el (ln| s |)GG* 
e -i§i e ie 2 _ e -is 2 ~t8s(F 1 ( a )-F 2 (u)) e -i(\n\s\)G'G 



Note that 
(5.4) 
Finally, we have 

s 



i 

s 

1 

J e 3*v+'^ 



^ (difi(a) (/3 1 ,e^l s D GG *a 1 ) w +a 2 / 1 (a)(e-(H S |) G * Ga2j/32 ^) 

l ai (ft/i(a) (e^I'D^'ai,^^ +ft/i(o) (/3 2 , e^'l^aa) 

+C(i? 2 s 



/ « 

2 „-2\ 



NONLINEAR LANDAU-ZENER FORMULA 



19 



Similarly, we obtain 

F 2 {u)u 2 = F 2 (a)e lf>2 a 2 + -F 2 (a)e lf)l (3 2 

s 



(d 1 f 2 (a) (e l ( ln l^ GG *a 1 ,/3 1 ) w +a 2 / 2 (a)(/3 2 ,e-^| S |)G- Ga2 ^) 
(dxMa) (/3i,e^l s l) G * G a 1 ) w +a 2 / 2 (a)(e-^| S |) G *G a2j/32 ^) 

+0{R 2 s- 2 ). 

By the definition of j3j, jj, kj and Kj for j G {1,2}, one gets that u(s) satisfies 

-id s u = V{s)u + F{u)u + 0{R 3 s- 2 ). 

Step 2. Existence of an exact solution with prescribed behavior. We first 
prove the result in the linear setting. Consider 



v-(s) 



e ts 2 / 2 +i(\n\s\)GG'/ 2ai 
e -ts 2 / 2 -i(\n\s\)G'G/ 2a2 



sq e R + such that R 3 s 1 < e, vq — us=o an d v E = (vf, v 2 ) solution to the linear 
equation 

-id s v £ = V{s)v e ; wf s= _ So = v-(s ). 

Then, by step 1, 

±\\v ( S )-v s ( S )\\ n *=O(R 3 s- 2 ). 
Therefore, there exists C > such that for s < —sq, 

\\v Q (s)-v E (s)\\ n 2 ^Ce. 

Consequently, 

(5.5) \\v-(8)-if(a)\\ n *^Ce. 

In particular, this implies that the scattering states (a 1 / 11 ' 6 , a 2 in ' £ ) associated with v e 
satisfy ||a^ m,e — aj\\u < Ce for j € {1, 2}. 

Let us now study the nonlinear case. We set 

J —oo 

This integral is well-defined, in view of (|5.3|) — (|5-4[) . since integration by parts yields 

$i(s) = O ( r-r ) as s -> -oo. 
\\s\J 

Define v by Vj = e~ l5sF ^ Q '~* 5 * J '^iij(s). In view of the first step, it satisfies 

-id s v = V{s)v + 0(R 3 s- 2 ). 

Moreover, by construction, 

\\v(s)-v-(s)\\ n 2 — > 0, 
s — y — oo 

so in view of 

±\\ v c(s)-v( S )\\ n2 =0(R 3 S - 2 ), 
as 

there exists a constant C > such that, for s < — sq, 

\\v(s)-v(s)\\ n 2^Ce. 
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Now set 

J — OO 

Then u e solves (|1.5j) . and since F is gauge invariant, 

^( s )=/ (F 3 (u%t))-F 3 (a))dt = o(J-^j as S ^-co. 

We infer 

(5.6) 3C>0, Vs<-s , \\u e {s) -u(s)\\ < Ce 

and the scattering states (af,a|) associated with u e satisfy \\a?j — ctjWn ^ Ce for 

./• {1,2}. 

Finally, if the unit ball of H. is compact, we can extract from u e a converging 
sequence which still solves (|1.5j) and whose scattering states are (cta,^)- 

6. Analysis of the Scattering Operator 



We now prove Proposition 11.71 and we take into account the dependence of the 
nonlinear term with respect to the parameter S. Because of (|4.9[) . we have cj lm = 
QrCj and 

u = cxp fiS (F(u(t)) - F(u)) dr^j w lin . 

Similarly, 

a = exp Us J (F(u(t)) - F(a j) dr^j a lin . 

In particular, we have F(u>) = F(w lln ) and F(a) = F(a hn ) since 
(6.1) || WJ || W = m in ||« and \\ aj \\ n = ||af \\ n for j E {1,2}. 

Besides, in view of 

u(s) = exp (iS f F{u{r))dT J u lin , 



we obtain F(u) — F(u hn ). We deduce 

u = exp (iS (F(u lin (r)) - F(uj lin )) dr^j cu l[n = c 

a = exp fiS J (F(u lin (r)) - F(a lin )) dr 



where A* = diag(A^,A^), and A^ are defined in (|1.11|) and (|1.12|) . respectively. 
Therefore, the nonlinear scattering operator S$ = S$(a) such that uj — Sga satisfies 

(6.2) S s (a) = e ^ A+ ( Q )^ lin ( e - i4A_ ( Q >a) , 

and the scattering operator Ss has the following expansion 

S s = S lin + iS (A+S lin - S lin A-) + O r {8 2 ) 

where O r (S 2 ) = E ie{1 , 2} (<W|A+|| 2 ) + 0(<5 2 ||A-|| 2 )). 
It remains to obtain an upper bound for A ± . We write 

rso r+oo 

A+ = (F j (u(s))-F j (u>))ds+ (F^u^-F^ds. 
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By the continuity of the functions f\ and fi (defined in (|1.6|l ) and because of 
II Willi/ + ||w2|||/ = llwiHl^ + ||w2|||/ = 1, we have 

\F j {u{s))-F i {^))ds ^fsup|/A so- 

Lemma |4~31 gives the existence of a constant C such that, 

H-oo 

(6.3) / (F.,M S ,z))-F>(z))Ks <Ci2V- 

■/so 

In order to optimize the upper bound, we choose sq = R whence 

O r {5 2 ) = 0(RS 2 ), 
which completes the proof of Theorem 11.71 

Appendix A. Rigorous derivation for a double- well potential 

A.l. Mathematical framework. We give more details concerning the derivation 
of p.5p in the case of a condensate in a double well. This is achieved by adapting 
the approach from [23) . from which several intermediary results are borrowed (see 
also [I] for some refinements to [23] in the confining, one-dimensional case). We 
shall derive the model (|A.3|) as an envelope equation in the semi-classical limit. 
Rewrite ()2.3j) in the presence of physical constants: 

(A.l) ih^ + -M h = V h (t,x)^ h + e(Ti)\^ h , 

where e(h) is a coupling constant whose value will be discussed later on. Note that 
we consider a slightly more general framework than in Section [2j x € R d , with 
d ^ 1. We assume that V h (t,x) — V s (x) + n(K)tV a {x), with V s and V a independent 
of h. 

We first describe the assumptions performed on the potential V s and discuss the 
first consequences. 

Assumption A.l (Symmetric potential). The potential V s € C°°(M ) is a smooth 
real-valued junction such that: 

(1) The potential V s is at most quadratic, 

d a V s e L°°(IR d ), V|a|>2. 

(2) V s is symmetric with respect to the first coordinate: 

V s (-x 1 ,x 2 , ...,x d ) = V s (xi,x 2 , ■ ■ -,x d ), Vx € K d . 

(3) V s admits two minima at x = x±, where x_ and x+ are distinct and sym- 
metric with respect to the first axis. Moreover, 

V s {x) > V s min = V(x±), Vx e R d , x ^ x±, 

and 

K min < liminf V s {x) =: V~. 

\x\— 5-00 

(4) The minima x + and X- are non- degenerate critical points: W(x±) = 
and V 2 V{x ± ) > 0. 

Remark A. 2. As noted in [23] . the last assumption (non-degeneracy) is not crucial. 
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We denote by 

h 2 

H =-—A + V s . 

The operator Ho admits a self-adjoint realization (still denoted by Hq) on L 2 (M. d ) 
(see e.g. [22]). Let u(Hq) = Od U cr css be the spectrum of the self-adjoint operator 
Hq, where ad denotes the discrete spectrum and tr ss denotes the essential spectrum. 
It follows that 

adCiVrW^) and = [V~, +00). 

Furthermore, the following two lemmas hold, which follow from [2]: 

Lemma A. 3 (Lemma 1 from [23] )■ For any h € (0, H*], for some h* fixed, it follows 
that: 

(i) ad is not empty and, in particular, it contains two eigenvalues at least. 

(ii) The lowest two eigenvalues A^. of Hq are non-degenerate, in particular, 
A^ < A_ . There exists C > 0, independent of h, such that 

X6a(/r )\[X»,A»] V ± ' 

Lemma A. 4 ([2 , and Lemma 2 from [23]). Let ip\ be the normalized eigenvectors 
associated to A±, then: 

(i) Lp\ can be chosen to be real-valued functions such that 

if% z (-xi,x 2 ,.-.,x d ) = ±tp%(xi,x 2 ,. . . ,x d )- 

(ii) (fl G EnL°°(R d ), where 

S = {/ e ff x (R d ), x h> \x\f(x) e L 2 (R d )}. 

(iii) There exists C independent of h such that for all h G (0, h*\, 

\\<A\\ lp(jl*) < Ch-i^), Vj»G [2,oo], 

||v^|| i2(R£i) ^ch- 1 ' 2 , ||^||| i2(Ri) < a 

The single- well states are then defined as in (I2.6[) . They satisfy: 
ip%(-xi,x 2 , ■ • • ,x d ) = y\(x\,x 2 , ■ ■ -,x d )- 

J I = O (e~ c / h } for all c <T, where F denotes the Agmon distance 



between the two wells 



F = inf / y/V s (x) - V™ ia dx > 0. 

7 connecting the two wells J 

• For any r > 0, there exists C > such that 

\<p h R (x)\ 2 dx = l + o(e- c / h ) , f \ V h L (x)\ 2 dx = l + o(e- c / h ). 

lB{x+,r) ^ ' JB(x-,r) ^ ' 

We denote by 

n* = i-«^,->^ + <^,->^) 

the projection onto the eigenspace orthogonal to the bi-dimensional space associated 
to A±. Finally, we define the splitting between the lowest two eigenvalues 

\ h \h 

(A.2) u h = — -y— L - 

Then for all c < T, uj h = O {e~ c / h ). 
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We now describe the assumptions made on the potential V a . 

Assumption A. 5 (Antisymmetric potential). The potential V a £ C°°(M ) is a 
smooth real-valued function such that: 

(1) The potential V a is bounded, as well as its derivatives, 

d a V a e L°°(R d ), V|a| ^ 0. 

(2) V a is antisymmetric with respect to the first coordinate: 

V a {-x 1 ,x 2 , ...,Xd) = -V a (xx,X2, ■ ■ -,x d ), Vie G K d . 

Remark A. 6. We note that Assumption IA.1I and I A. 51 are satisfied in the case (|2.5[) 
borrowed from [17] . 

A. 2. An approximation result. In this section, we prove: 
Proposition A. 7. Let d ^ 2. Let V h be as in Section \A.l[ with 

(uj h ) 2 

k = k(H) = V—p—, 
for some rj € R independent of K. Suppose that e(H) is given by 

e(fi) = 8ui h h d / 2 , 

where 5^0 does not depend on h. Suppose also that the initial datum is of the 
form 

ip h (0, x) = a[ J ip r l(x) + a R ip R (x), ctL, otR G C independent of h. 
Define the approximation solution V'app by 

(u>^t\ / 0J^t\ 

— ) <f%(x)+a R i—j <pk(x), 

where (aL,a R ) = (ol(t), <2.r(t)) solves 

{id T aL = r]TOL - an + 5 h \aL\ 2 aL] a L \ T=0 = oll, 
id T a R = -r)Ta R - a L + 5 n \a R \ 2 a R ; a fl | T=0 = a R , 



5 h = 5h d ' 2 I <p\ = 8h d ' 2 I <p R 



(A.3) 
and 



is uniformly bounded in h £ (Q,h*]. Then we have the following error estimate: 
there exist c, C independent of h such that for all 7 < T, 

sup U h (t) - e-*< A -+ A +>/( 2ft >< p (i)|| ia < C^l\ 

\t\^es/h /u h 

Remark A. 8. The case d = 1 and S < could be considered as well, leading to the 
same result. Considering this case would just make the proof a bit longer, and we 
refer to [23] for the adaptation. 

Remark A. 9. Since the range for the slow time r = ui h t/h is —c/Vh ^ r ^ cj\fh, 
the above approximation result is a large time result, which is consistent with the 
large time study of (]A.3|) . or, more generally, of (jl.5l) to which one reduces thanks 
to the change of variables s = ^rjr (we then have G = —r)~ x l 2 and SF(u) — 
7y _1 / 2 5 ft Diag(|tti| 2 , |u2| 2 ))- When h is small, Theorem 11.31 gives asymptotics for the 
profiles ql and a R of ^a P p f° r large times < cVh/uj h ) and the connexion between 
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the profiles for t < and t > involves the Landau-Zener transition coefficient e 
(at leading order when h goes to 0). 

Proof. For simplicity, we write the proof for t and it naturally extends to t 0. 
Step 1: Preliminaries. We begin by proving estimates on ij; h , then we perform 
the rescaling suggested by the form of the approximation solution. First, it follows 
from Lemma TA.4I and [4] that for fixed h > 0, (|A.1[) has a unique, global solution 
ip h g C(K; S) (S is defined in Lemma [A.4|) . In addition, if we set 

Mass: M h {t) = U h {t)f L „ 
Energy: E h (t) = \\\tN ^ h (t)\\h + ^W^WWh + / V(t,x)\^ h (t,x)\ 2 dx, 

then we have the a priori estimates: 
dM h 

u . 

d t V(t,x)\ip h (t,x)\ 2 dx = n(h) { V a (x)\^ h (t,x)\ 2 dx, 



dt 
dE h 
dt 



where we have used (|2.4I) for the last equality. Note that (|A.1[) is not a Hamiltonian 
equation, unlike the one studied in pQ, where the Hamiltonian structure is crucial. 
The conservation of mass and the form of the initial data yield 

M\t) = U h {t)\\ L , = ||V/ l (0)|| L2 = M h (0) < C, 

for some C independent of h. We have 

E h (0) = {a LV h L + a R v% a L H <p h L + a R H <p%) + ^-\\a L <p h L + a a <p%\\i<. 

Introduce fl h = (A* + \\)/2 and notice that fl h = V™ in + 0{h) by Lemma[0] 
Noting the identities 

(A.4) H ^=nVE-w% Ho^ = ^-<A>L 

we infer from Lemma IA.3I and Lemma IA.4I that 

E n (0) = V s min M H (0) + O(h). 
Therefore, since V a is bounded, 

|| W n (i)||| 2 s? 2E h {t) - 2V s min M h (t) + Cn{h)tM h {t) 

^2E h (0)-2V> riin M R {Q) + 2 K {h) f f V a (x)\ip h (t, x)\ 2 dx + Cn{h)t 

JO J R d 

(u h ) 2 

h 

Next, as in [53], we consider the slow time 

T = — > 0, 

n 

and the new unknown function 

y h (r, x) = ^> h (t, x y nH i h = ip n (t, x y nhT ^ h . 

It solves 

(A.5) id T V H = i (H a - fl h ) f fi + r]TV a ^ h + 5h d ' 2 \^ h \ 2 y h , 
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and in view of the above estimates, 

= n(~\\- llV7\[r ft ^ll?„ < n I — u 



(A.6) ||**(r)IU' - ||* ri (0)|| L2 = 0(1); || W\r)||i 2 < C ( \ + 



We decompose SSf h as 

(A.7) * R = / l + ^, - n*¥ n , 

so we can write </? as 

V \T,x)=a.l(T)<d(x) + *&T)<p%(x), 

for some complex- valued coefficients a| and a^. Projecting (IA.5I) . and using (IA.4I) . 
we find: 

ikl = -z n R + r)T [ v a ^l + sn d/2 [ |* B | 2 *Vi, 
ik% = - a n L +rrr [ K*Vh + ^ /2 / l* B | a *V«, 

The proof of Proposition I A. 71 consists in showing that , R are close to a^/n, and 
that is small. This is achieved in two more steps. 

Step 2. A priori estimates on a|,, a^ and ip*. By definition, we have a^ = 
L<j ty h y>\, so Cauchy-Schwarz inequality yields 

ll/(r)|| 2 L2 = |ai(r)| 2 + |a* (r)| 2 < (M h f (||^|| 2 L2 + ||^||| 2 ) < 2 (M*) 2 . 
Decompose the nonlinearity acting on ^> h as 

\-q> h \ 2 -q, h = \(p h \ 2 (p h + R H . 
Using the a priori estimates on a. h and Lemma IA.4| we have 

IH^iVL < ii^iii-ii^iu. < c (ii^nioo + n^ni«) < ch- d '\ 

Since we have the pointwise estimate 

|R"KC(|^| 2 |^| + |^| 3 ), 

we infer 

ll^lU^C^-^H^II^ + ll^llla), 

and Gagliardo-Nirenberg inequality yields (d ^ 2) 

IIV'c lli 6 (R d ) ^ CIIV'c llx,2(ffid)IIWc lll 2 (IR d ) ■ 

Since d ^ 2, we can factor out HV'c lli 2 (R ti )' an d (|A.6p gives 

u n c \\l* m < cn^ii L2(R ,)fi- d/2 ^i + (^p) 1 ) . 

Therefore, 

(A.8) ||R R || L2 < Ch- d ' 2 (l + (^) ' ) ||^|| £3 , 
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CO T X 



and IvI/"! 2 ^" satisfies a similar estimate. We infer 

|a£(r)|+|afc(r)| ^ C (l + r) + C (l + {j- 

Since d ^ 2 and uo h = C(e- c / fi ), we can simplify the above estimate: 
14(01 + |4M| < C(l + r), Vr^O. 

From this we infer 

(A.9) ||^1 i2 < |4(r)| U^IU* + |a|(r)| ||^|| L2 < C(l + r), 

and, again from Lemma lA.41 

(A.10) ||9 r (|^|V)L 2 < n<P*\\l- \\d T ^\\ L2 ^Ch- d / 2 (l + r). 

Step 3. Stability estimates. In view of (|A.7j) . we first prove that ip^ is small. 
Since ijA T _ Q = 0, Duhamel's formula yields 

(A.ll) $( T ) = -if] f T e -*(H°-n R )(r- S )/^ ( gI1 h (y^Hj ds 

Jo 

- i8h d ' 2 f T e -*lH -n*)(T-s)/^ n h (|$»|2 $ »j (s)ds , 
Jo 

Each term is treated in a similar fashion: when \E fft is replaced by y> , we perform 
an integration by parts, and for the remaining term, we use directly the a priori 
estimates. For the first part of (|A.11[) . we write 

rr} (v a * H ) = rr} (v aV h ) + nj (v B v£) , 

and set 

Jo 

4 l (r) = f T c -*(H -n»)(r-s)/^ ( sIl h (y a ^) ( s )) rf s . 
Jo 

Integrating by parts, 

iHr) = ^ e -i(H«-^)(T-s)/^ , Hq _ n «N-i n fi / y »n (s) T 

o 

Jo 

From Lemma [A. 3 [ there exists C independent of h € (0, ft*] such that 



^ c. 

L 2 ^L 2 



We infer, using (|A.9|) . 
For , we have directly 



For the nonlinear term in Duhamel's formula (the second term of (|A.11|) ), we also 
write 



\%(t)\<C s\\^(s)\\ L *ds. 
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and set 



If (r) = %m / e -i(Ho-n ,l )(r- a )/«» n » (|^| V) ( s ) ds , 
Jo 

i!(r) = h d ' 2 r e -(a-si s )(- s )/- fi n c ¥( s )d s . 



\^( S )\\ L 2dS, 



We have, directly from (|A.8|) . 



and performing an integration by parts for 1$ , using (|A.10[) . we have 



|/f(r)|<C^-(l + r 2 ) 



Since d ^ 2 and a; decays exponentially, we come up with: 

ll^(r)|| L2 < C (j- (1 + r 3 ) + jf (1 + S )||^( S )|| i2 d S ) 



Gronwall lemma yields 

ll^(r)|| i2 <C^(l+r 3 )e c ^ 2 ). 

Recalling again that uj h decays exponentially, we can write that for all Co < T (the 
Agmon distance between the two wells), 

This is small as h — > 0, provided that t 2 <C /h. More precisely, there exist c\, C2 > 
independent of ft such that 



(A.12) 



To conclude the proof of the proposition, set 



Subtracting the equation for from the equation for a^, we have 

JR d 

sn^ [ (|^| 2 ^-| aL |V|^iM)^. 



We have 



/ V a (* h -a L v n L )vl= f K(a|^ + ^+w^|)^, 
therefore, since the product V'iVfl. decays exponentially in ft, 

K - aL <p*) <p* < C (c- c / R + ||V4|| L2 + |w||) . 
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A similar estimate can be established for the other source term in the equation for 
w^. The equation for is handled in the same fashion, and using (|A.12|) . we end 
up with: 

|w£(r)| + |4(r)| < C (|v*(r)| + |v* (r)| + e~ c ' n ) , < r < 
Gronwall lemma and (|A.12[) then yield Proposition I A. 71 □ 
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